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Oh, Abstract. In the present article magnetic Laplacians on a graph are analyzed. We provide a 

. complete description of the set of all operators which can be obtained from a given self-adjoint 

Laplacian by perturbing it by magnetic fields. In particular, it is shown that generically this set 
is isomorphic to a torus. We also describe the conditions under which the operator is unambigu- 
ously (up to unitary equivalence) defined by prescribing the magnetic fluxes through all loops of 
^ ' the graph. 

^ ! 1. Introduction and Main Results 

O . 

O ■ Magnetic Laplacians on finite graphs appear in a number of physical applications. The major 

interest in operators of this type originates from the study of quantum transport in mesoscopic 
ps| \ networks (see, e.g., [1], [3], [4], [5], [6], [7], [12], [15]). Magnetic Laplacians have also been 

O ■ the subject of several studies in the context of quantum chaos [24], [25], [26]. Graph-theoretical 

generalizations of the Harper operator provide discrete models describing the dynamics of a 
Oh' quantum particle in the presence of a magnetic field (see [34], [28] and references therein). 

. The most intriguing feature of Laplacians on graphs is the relation between their spectral 

properties and the geometry of the graph. For discrete Laplacians the study of this relation is 
one of the central issues of spectral graph theory (see, e.g., [16]). As an example, the multiplicity 
of the zeroth eigenvalue of the discrete Laplacian equals the number of connected components 
of the graph. In the presence of an external magnetic field the effect of the topology of the graph 
^ ■ becomes more pronounced. Thus, the quantum conductance in networks is known to be related 



^ ; to Chem numbers [3], [4], [5], [6]. 



In the present article we study general differential self-adjoint magnetic Laplacians on finite 
graphs. This work is a continuation of our previous studies [20], [21], [22]. 

We consider an arbitrary nontrivial connected (metrical) graph S with a finite number « > 
of external lines and a finite number ni > of internal lines (henceforth also called edges), 
m + n^Q. More precisely, this means that outside of a finite domain the graph is isomorphic to 
the union of n positive half-lines. Any internal line ends at two, not necessary different vertices 
and has a finite length. We assume that any vertex z; of S has non-zero degree deg(z;), i.e., for 
any vertex there is at least one edge (internal or external) with which it is incident. 

Let the set £ (j£j = n) label the external and the set J = m) the internal lines of the 
graph. We assume that the sets £ and J are ordered in an arbitrary but fixed way. To each 
e G £ we associate the semi-infinite interval [0,oo) and to each / G J the finite directed interval 
[0,a,], where a, > is the length of this line. With this association the graph becomes directed, 
such that the initial vertex of an edge of length a, corresponds to x = and the terminal vertex 
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corresponds to x = a,-. The external lines are assumed to be directed in the positive direction of 
half-hnes. We denote by a the collection of the lengths {a,},£j. 

Though not necessary but rather as a motivation for our approach we temporarily assume 
that the graph is isometrically imbedded in three dimensional Euclidean space M?. Given a vec- 
tor potential A G (C(M^))'', whose curl gives a magnetic field, the question arises what is the 
reasonable self-adjoint operator describing the evolution of quantum states on the graph in the 
presence of this external magnetic field. A common intuitive construction goes as follows. We 
assume that all external hues and all edges of the imbedded graph S are oriented smooth curves. 
Let Tj{x) (y G £ U J) denote the unit tangential field on the j-th edge of the graph chosen in ac- 
cordance with the orientation of the edge. We set Aj{x) = {Zj{x),A{x)) where (•, ■) denotes the 

inner product in M^. The vector potential leads to the differential expressions (— i^ — yiy(x))^ 
on the external lines or edges of the graph with appropriate boundary conditions at the vertices. 
More precisely, consider the family i/a = {Yjjjeeud of complex valued functions defined on 
[0,00) if y G £ and on [0,a,) if j G J. Formally the Laplace opertor is defined as 

(1.1) {A(A,B;A)xif)jix)=( — -iAj{x)j Yji^), jeJuE 
with the boundary conditions 

(1.2) AY + Bi^-iAY)=0. 
Here 

\{Viiai)}ieoJ \{-¥liai)}ieO^ 
and A is the diagonal matrix 

(1.3) A = diag{{Ae{0)}ee£,{M0)}ie3,{-Mai)}ie3) 

with the same ordering as used for i/a and i/. If yi = we will simply write A(A,S) instead of 
A(A,B;0). ~ ~ 

Actually we do not need to imbed the graph in M? and we may simply prescribe the vector 
potentials A{x) = {Aj{x)}j^j(j£^ on all external Unes and edges of the graph. So having given 
this motivation, our further discussion will be completely intrinsic. 

Below we will prove (Theorem 2.1) that the operator A{A,B;A) is self-adjoint if and only 
if AB^ is symmetric (with t denoting the Hermitian conjugation) and the {n + 2m) x (n + 2m) 
matrix (A,S) has a maximal rank. This is an extension of a result in [20] to the case of magnetic 
Laplacians. 

The local gauge transformation 

(1.4) G: Yj{x)^\ifj{x)exp{iXj{x)}, jeJu£ 
with 

(1.5) Xj{^)=Aj{x) 

eliminates the vector potential in the differential expression (1.1) and in the boundary conditions 
(1.2). However, the gauge transformation changes the boundary condition: It transforms the 
magnetic Laplacian A{A,B;A) to the Laplacian A{AU,BU) without magnetic field, 

G-^A{A,B;A)G = A{AU,BU), 
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where the transformation U = Kg is given by the diagonal, unitary matrix acting on the boundary 
values Y and y' 

(1.6) U = diag({e*(»l}^^^,{e*<»)}_.^^,{e*('")}J 
with 

XM) := Xi{0) + r Ai{t)dt, i G 3. 
Jo 

In the other words, the perturbation of the Laplacian A{A,B) by the vector potential A is 
equivalent to the transformation of the boundary conditions {A,B) {AU,BU) by some diago- 
nal unitary matrix li. Conversely, given a diagonal unitary {n + 2m) x (n + 2m) matrix It there 
is a vector potential A (possibly zero) such that A(A'Lt,B'Lt) can be obtained from A{A,B;A) 
by means of a local gauge transformation. We will call A{AU,BU) a magnetic perturbation of 
A{A,B). 

This correspondence can be described more precisely in a group-theoretical setting. Let 
U = U(9) denote the group of all diagonal unitary {n + 2m) x [n + 2m) matrices enumerated 
by the elements of the sets £ and 3 in the same ordering as in (1.6). Let Uq be its subgroup 
consisting of those matrices with the structure (1.6) with exp{i;|^, (0)} = exp{i;|^,(a;)} for any 
i G J. Equivalently, the elements in Uq correspond to gauge transformations (1.4) with vanishing 
vector potential A. With this notation the cosets U /Uq of Uq in U are, obviously, in a one-to-one 
correspondence with the points of the m-dimensional torus T'" given as 

?, = exp|iy^ G S, j G J. 

The first main result of the present work provides an answer to the following question 

• What is the orbit O (A, S) := {A(AU,S1()| U G U} of the operator A(A,B) in the set of all 
self-adjoint Laplacians under the action of the group U? Roughly speaking, the question 
is how many different magnetic perturbations of a given Laplacian A(A,6) there are. 

The orbit of A(A,B) under the action of U is isomorphic to the factor group U/W(A,S), 
where W(A,5) is the isotropy group (or stabilizer) of the point A{A,B), i.e. the subgroup which 
leaves A(A,B) invariant. Recall that by a result of [20] A(A1X,S1X) = A(A,S) if and only if there 
is an invertible {n + 2m) x (« + 2m) matrix C such that CA = Alt and CB = BU. 

The following theorem provides an explicit description of W(A,S). Let A* denote the Moore- 
Penrose pseudoinverse of A (see, e.g., [20] for a short presentation of the main facts related to 
this notion). 

Theorem 1. The element It G U belongs to W(A,B) if and only if 

(i) It leaves both KerA and KerB invariant; 

and 

(ii) one of the following equivalent commutator relations 

(1.7) [U,A*B]z = for any z G fi-^RanA nRanS) 
or 

(1.8) [U,B*A]z = for any z G A-i(RanAnRanB) 
isfiilfilled. 

In particular, j/RanA n RanS = {0} then A(A'Lt,B'Lt) = A(A,B) for allUeU such that U 
leaves both KerA and KerB invariant. 
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Note that since U is unitary, the condition (i) implies that both orthogonal decompositions 
KerA [KerA]-^ and KerS [KerSj-^ of 0"+^'" reduce U. Also we remark that KerA n KerS = 
{0} since {A,B) has maximal rank. 

There are examples of boundary conditions (A,S) (see, e.g.. Example 2.7 below) where 
W(A,B) = U and, thus, the orbit 0{A,B) consists of the only point A(A,S). 

An important class of boundary conditions are local boundary conditions. They are such that 
they couple only those boundary values of Y and of its derivative ijf' which belong to the same 
vertex. The precise definition is as follows. 

Let S{v) C £ U J denote the star graph of the vertex v eV, i.e., the set of the edges adjacent to 
V. Also, by 5_(i;) (respectively ^^(i;)) we denote the set of the edges for which v is the initial 
vertex (respectively terminal vertex). Obviously, S+{v) nS-{v) = if S does not contain a 
cycle of length 1, i.e. a tadpole. 

Assume the elements z of the linear space C""*"^"" are written as 

(1.9) Z={{Ze}ee£,{^~^}ieO,{i'^^}ie3f- 
Consider the orthogonal decomposition 

vev 

with being the linear subspace of dimension deg(t;) spanned by those elements (1.9) of 

,^n+2m ^IjI^Ij gjjjigfy 

Ze = for all ee8.\S{v), 

(1.10) 4~^=0 forall ie3\S-{v), 

zS+^=0 for all ie3\S+{v). 

Definition 2. Given the graph S, the boundary conditions (A,S) are called local if there is an 
invertible [n + 2m) x (« + 2m) matrix C and linear transformations Ax, and Bd in such that 
the simultaneous direct sum decompositions 

(1.11) CA = 0Ap and CB=^B^ 

vev vev 

are valid. 

Otherwise the boundary conditions are non-local. 

As outlined in [20] (see also [27]) for an arbitrary self-adjoint boundary condition (A,B) there 
is always a graph with external lines labeled by £, internal hnes labeled by J with a set of lengths 
a, for which these boundary conditions are local. This graph is unique under the requirement 
that the number of vertices be maximal. We will elaborate on this in [23]. 

Consider the group Wq consisting of those elements U of U which admit a decomposition 

XL = Uj), Ud : — > fiiv 

vev 

with U^; being a multiple of the deg(i') x &eg{v) unit matrix (see Definition 2.8 below). It is 
an immediate consequence of Theorem 1 that Wq is a subgroup of the isotropy group W(A,S) 
for any local (self-adjoint) boundary conditions. For non-local boundary conditions it may well 
happen that W(A,S) C Wq (see Example 3.4 below). 
The next question we turn to is 

• What is the correspondence between the points in the orbit of the operator A(A,S) under 
the action of the group U and the magnetic fluxes through the loops (i.e., closed paths) 
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of the graph S? Given a Laplacian A{A,B) can its magnetic perturbation A{AU,BU) 
be unambiguously defined by prescribing the magnetic fluxes through all loops of the 
graph? 

We will give an affirmative answer in Theorem 3 below. To formulate this result we need 
some additional notation. 

A vector potential A = {Ajjyegua defines a flux map, by which we associate to any loop y in 
the graph S the complex exponential of the magnetic flux through 7, 

(1.12) ^Air) = exp i^i(j^sign{i) £' Mx)dx^ | . 

Here the sum is taken over all edges / G 3 constituting the loop 7 with sign(/) = +1 if the 
orientation of 7 coincides with that of the edge / and sign(/) = — 1 otherwise. If the graph S is 
imbedded in M?, by Stokes theorem the expression in the brackets (...) has indeed the meaning 
of a flux through any oriented bordered manifold M with boundary dM = 7. 

By established equivalence between vector potentials and transformation matrices U G U the 
flux map (1.12) can be alternatively expressed in terms of the matrix elements of IX (1.6) as 

(1.13) ^u(r) = exp i^i{Zsign{i){Xi{ai)-Xim) | 

and where the sum is as in (1.12). 

Loops of the graph S can be realized as elements of the additive AbeUan group Hi{S,^), 
the first homology group of the graph S. Recall that this group is the Abelianization of the 
fundamental group Its elements are formal Unear combinations with integer coefficients 

of the "basis" loops in the graph, 

p 

For any U G U the flux map <I>xx defined by (1.13) can uniquely be extended to a group 
homomorphism from //] (S, Z.) to S. The set of all such homomorphisms forms an Abelian group 
denoted by Hom(//i(S,Z),S). The map p : IX 1— > is obviously a group homomorphism. 
Actually, it is an epimorphism (Lemma 3.2). It is clear that Uq ^ Ker p. In Lemma 3.3 below 
we will prove that Ker p = UoWq. 

The second main result of the present work is given by the following theorem. 

Theorem 3. Assume that the self-adjoint boundary conditions {A,B) are local. Then the groups 
\J/{Uo\N{A,B)) and}iom{Hi{S,Z),S)/p{\N(A,B)) are isomorphic. In particular, if\N(A,B) = 
Wo, then the groups U/(UoWo) an<i Hom(77i(S,Z),S) are isomorphic. 

Whenever the local boundary conditions {A,B) are such that \N(A,B) = Wq, Theorem 3 states 
that for any given magnetic flux there is a unique set (up to transformations defined by elements 
of Uo and Wq) of phase factors IX G U giving this flux. The situation changes drastically for 
boundary conditions such that W(A,S) is strictly larger than Wq. The transformations (A,B) 1— > 
{AUjBU) with elements IX of W which are not elements of Wq do not change the operator 
A{A,B). However, they change the magnetic flux! In the latter case, however, there is still a 
subgroup of Hom(//i(S,Z),§) corresponding to those fluxes which remain unaffected by all 
transformations {A,B) i—^ (Alt, BIX) induced by the group W. This subgroup may be realized as 
a factor group Hom(//i (S, Z) , S)/p (W). Example 2.7 below shows that W may be as large as 
U such that Hom(7/i(g,Z),S)/p(W) is ti-ivial in this case. 

According to Theorem 3 magnetic Laplacians depend on the magnetic flux through any loop 
of the graph 2;r-periodically. In the special case of a ring (in our context a graph formed by a 
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single internal line with coinciding initial and terminal vertices) this fact is well known in the 
physics literature as the Byers-Yang-Bloch Theorem [9], [8]. 

The following corollary of Theorem 3 states that if S is a tree then an arbitrary magnetic 
field A{x) does not change the Laplacian A(A,S) for any local boundary conditions (A,B) in the 
sense that there is a local gauge transformation G such that ^{AILg^BUg) = A(A,B) with Ug 
being given by (1.6). The precise statement is as follows. 

Corollary 4. Tjf ;ri(S) is trivial (and hence /?i(S;Z) is also trivial), i.e. if the graph S contains 
no loops, then U = \JoW{A,B) for arbitrary local boundary conditions {A,B). 

The requirement of locality is crucial since non-local boundary conditions may have the same 
effect as a loop. An example of such situation is presented in Example 3.4 below. 

As already mentioned magnetic Laplacians on graphs appear as models of different physical 
systems. In particular, the Chalker-Coddington network model [15] was designed to describe the 
semiclassical motion of a single electron in the presence of a perpendicular uniform magnetic 
field and a random potential. Thus, in this model the internal lines of the network (graph) are 
determined by the equipotential Unes of the potential and the vertices by its saddle points where 
two equipotential lines closely approach one another. Correspondingly, the tunneling and trans- 
mission probabilities then give rise to boundary conditions at the vertices. As a consequence of 
the results of the present work a full quantum version of the Chalker-Coddington model can be 
developed. 

Acknowledgement. One of the authors (R.S.) would like to thank M. Schmidt and E. Vogt for 
advise and help. Discussions with H. Schulz-Baldes have been very helpful. 

2. Laplacians With Magnetic Field 

To the triple (£,5, a) with £, J, and a as being given in the Introduction we associate the 
Hilbert space "K = 5{(£, J,a) given as the orthogonal sum 

where "Ke = L^(0,oo) and J{/ = L^(0,a,). The iimer product in is given by 
{<I>,¥):k=Y. (t>e{x)xi/e{x)dx + y ^i{x)\ifi{x) dx. 

Elements of "K are written as column vectors 

(2.1) XI, = (^^^^^^^ , Wee^Ke, Wi&^i- 

Similarly we define the Sobolev space as 

W2'' = 0^2-2(0, oo) © 0w2'2(O,a,), 

where W^-^{Q,oo) and W^-^{Q,ai) are the usual Sobolev spaces of square integrable functions 
whose distributional second derivatives are also square integrable (see, e.g., [32]). We observe 
that the Hilbert space is independent of the particular graph constructed out of the set of data 
(£,J,a). 

Theorem 2.1. Let A{x) = {Aj{x)}j^'j\j£_ be continuous and bounded. Any two (n + 2m) y.(n + 
2m) complex matrices A and B satisfying 

(2.2) AB+-SA+ = 
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and such that the (n + 2m) x 2{n + 2m) matrix (A, 5) has maximal rank equal to n + 2m, define 
a self-adjoint magnetic Laplacian 

{A{A,B;A)y)j (x) =(J^- i-Ay w) Vji^), jeJuE. 

in IK corresponding to the boundary condition 

(2.3) Axi/+B{xj/ - iAxj/) = 

and with the diagonal (n + 2m) x (« + 2m) matrix A given by (1.3). 

Conversely, any self-adjoint magnetic Laplacian corresponds to the boundary condition (2.3) 
with some matrices A and B satisfying (2.2). 

Proof. Consider the symmetric operator defined by 



(A».(x) = (^^-iyi,(x)) Wjix) 



with domain 2)(A^) C W^'^ consisting of functions which vanish at the vertices together with 
their first derivative. It is clear that A^ has defect indices {k,k) with k = |£| +2|J|. 
On W^'^ we consider the Hermitian symplectic form 

with Aji being considered as a formal differential expression 

Let[U:W 
[y]a given as 



Let [ ]^ : ^^'2 c'2-{n+2m) jjjg surjective linear map which associates to each y the element 



A 



\ 



{Ve{0)}eee 

{¥imie3 
{Y'e{0)-^MO)¥eme€e 

{V//(0)-iA(0)v^(0)},ea 

V{- v^(«/) + ^Mi) wM)}i& J 

If A{x) is continuously differentiable then by means of partial integration one verifies 



v' - ^Av J 



^a{(I>,W) = L U0)¥'e{0)-(l>mWe{0)-2LAemMxifei0) 



(2.4) 



+ L - (l>mxifii0)-2iAi{0M0)xifii0) 



- L ^ii'^iWM) - {adWiiad - 2iA-(a/)0,(a,-) v/,-(a,-) 
={[(!>] A, Jiw] a) c2('>+2'") ) 

where 

■ I' 

is the canonical symplectic matrix on C'^i"+^>")_ Here I is the (« + 2m) x {n + 2m) unit matrix. 
In the general case we again obtain (2.4) if we approximate A{x) by continuously differentiable 
functions. 

Let the linear subspace 3VC(A,S) of C'^i"+^'") be given as the set of all 



A 



1 

xif' -iAxir 



j^2{n+2m) 
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satisfying 

AvA+B(v/-i^V^) = 0. 

All self-adjoint extensions of AJ^ are described by maximal isotropic subspaces M{A,B) [20]. 
By Lemma 2.2 of [20] we obtain the claim of the theorem. □ 

Remark 2.2. The technique of using Hermitian symplectic forms in extension theory seems to 
have been well known for a long time. The earliest reference we are aware of is [19, Section 
10]. In the context of differential and difference operators on graphs and in similar contexts this 
technique has also been used in [31], [30], [18], [11], [29]. 

We also mention the articles [13] and [14] where the self-adjointness conditions (2.2) were 
proven (without using the formalism of Hermitian symplectic forms) for the cases m = 0, n 
arbitrary and m = I, n = 0, respectively. 

The assumption that the vector potential is bounded formally excludes the case of homoge- 
neous magnetic field. We are allowed, however, to consider magnetic fields which are constant 
on an arbitrarily large finite domain. The behavior of the vector potential on the external Unes 
of the graph away from the vertices does not influence the operator A{AU,BU) obtained from 
A{A,B,A)hy a gauge transformation. Therefore, constant magnetic fields can be considered as 
well. 

Remark 2.3. We may replace the pair (A,S) by the pair {CA,CB) where C is any invertible 
(n -|- 2m) X {n + 2m) matrix. In fact such a replacement does not change the linear conditions 
(2.2) and (2.3) is still satisfied. In other words, the maximal isotropic subspace is left unchanged 
and so we have A{CA,CB;A) = A{A,B;A). 

Let x{x) denote an arbitrary family {Xj{x)}je3iiE of continuously differentiable real valued 
functions on [0,aj] if j G J and on [0,o°) if j G £. To such X we associate a unitary map G{x) 
in J{by 

(G(;c)VA);(x)=e'^^-Wv^X^)> ;e5u£. 

We claim that 

G{xYa{A,B;A)G{x)=A{AU,BU;A-x') 

with 

dxjixl 
dx 

and It defined by (1.6). Indeed, setting y = U<p we obtain 

AY + B{\/ + iAY) =AU^ + BU{^-i{A-x^)(j>)=0, 
where x' is the diagonal matrix given by 

/ = diag {{ximeeE, {Xm}i&, {-X- («/)},€a) 
and A is defined by (1.3). Choosing = A we obtain 

G{xyAiA,B;A)Gix) = A{AU,BU). 

We turn to the proof of Theorem 1 . For the proof we need the following auxiliary result. Re- 
call that A* denotes the Moore-Penrose pseudoinverse of A. Let Px be the orthogonal projection 
onto the linear subspace X C C""'"^'". Since the matrix (A,B) has a maximal rank, PRanA + PRanB 
is invertible. 
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Lemma 2.4. Assume that It leaves both KerA and Ker5 invariant. Then there is a matrix C 
satisfying 

(2.5) AUA* = CPr^a and BUB* = CPr^b, 
if and only if 

(2.6) AViA*y = BUB*y 

holds for all y G RanA PI RanS. IfC exists, it is invertible and given by 

(2.7) C = {AUA*+BUB*){Pr^a + Pr^b)-' 

Proof Assume there is a matrix C satisfying (2.5). Then (2.6) and (2.7) are obvious. We claim 
that C is invertible. To see this observe that from (2.5) and the fact that U leaves KerA and KerB 
invariant, it follows that C maps RanA and RanS onto themselves. Thus, the range of C is all of 

(^n+2m 

Now assume that (2.6) holds for all y G RanA fl RanS. We prove that the matrix C given by 

(2.7) satisfies (2.5). Since KerA* = KerA+ we have A* =A*PRsnA- Similarly, S* = S*/^anB. 

Therefore, 

(2.8) CPRanA = AUA*PRanA(^RanA +^RanB) ^ PrwA + BUB^PR^^aBiPRmA + Prioib)'^ PRanA- 

By a result of [2], 

PkanAiPRanA + Prwib) ^PkanB = PRanB{PRanA + PkanB) ^PriuiA = 2^anAnRanB- 

From this we immediately obtain 

^anA(-PRanA +-f'RanB) '-fkanA = PRanA — ^-^RanAnRanS- 

Thus, from (2.8) it follows that 

CPrwA = AIXA* — -AUA^/iianAnRanB + ^^^^^-^anAnRanB 

= AUA* - ^Al[A*/^anAnRanB + ^AUA*/^anAnRanB = AUA*, 

where we have used (2.6). The second relation in (2.5) is proved in the same way. □ 

Proof of Theorem 1. By Remark 2.3 A{AU,BU) = A{A,B) if and only if there is an invertible 
(n + 2m) X (n + 2m) matrix such that 

(2.9) AU = CA and BU = CB. 

Assume that these relations hold with some C. Then the first relation in (2.9) implies that 
Ally = for any y G KerA. Thus U leaves KerA invariant. Similarly, by the second relation in 

(2.9) KerB is left invariant by U. 

Since RanA* = RanA^ any y G [KerA]^ can be represented as j = A*z with some z G C"^^'". 
Therefore, from (2.9) it follows that 

(2.10) AUA* = CPR^nA and BUB* = CPr^b- 
Using Lemma 2.4 we conclude 

(2.11) AUA*y = BUB*y 

for all y G RanA nRanS. Multiplying (2.11) by A* from the left and using A* A = PRanAt we 
obtain 

(2.12) P^AtUA*y=A*BUB*y for any y G RanA nRanS. 
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Since RanA* = RanA' and since U leaves RanA' invariant, we can omit the projection on the 
l.h.s. of this equation. 

Let z = B*y. Since KerS* = KerS^ we have zj^O. It easy to see that y = Bz. Inserting this 
into (2.12) we obtain 

UA*Bz =A*BUPr^b^z =A*BUz 

for all z e S*(RanA nRanS) = B-i(RanA nRanS). This proves (1.7). The proof of (1.8) is 
similar. 

Conversely, let (i) and (ii) be valid. By Lemma 2.4 the matrix C defined by (2.7) satisfies 
(2.9). □ 

For given boundary conditions (A,S) let W(A,S) be the set of all 11 G U such that A{AU,BU) = 
A{A,B). This set, by construction the stabilizer of the point A(A,B), is obviously a group by 
well known arguments from Group Theory. However, we prefer to give a direct proof. 

Proposition 2.5. W(A,B) is a group. 

Proof. Assume that l(i,l(2 G W = W(A,S). By Theorem 1 U = U2U1 leaves both KerA and 
KerS invariant. In addition this theorem imphes that 

(2. 1 3) UiA*Bz = A*BUiz, UiA^Bz = A*BU2Z 

is vaUd for all z G B~^(RanA flKerS). Multiplying the first equation by U2 from the left we 
obtain 

UalXiA^Sz = U2A*BUiz. 

We claim that IXiz G B~^(RanA PlRanB). Then we will obtain from the second equation in 
(2.13) the equality U2UiA*Bz=A*BU2Uiz which implies that U2U1 G W. 

To prove that Uiz G S~^(RanA flRanS) it is sufficient to show that BUiz G RanA. Assume 
there is z G (RanA n RanB) such that BUiz -L RanA. Then from the first equation in (2.13) 
it follows that A*Bz = 0, i.e., Bz -L RanA, a contradiction. 

The proof that U^^ G W is similar and will therefore be omitted. □ 



For local boundary conditions the group W can be factorized as a direct product of its sub- 
groups Wt,, 



W(A,B)= X \NviAr;,Bv). 

veV 



(2.14) 

Example 2.6. An important example of local boundary conditions is given by the matrices 



(2.15) 



At, = 



/I 


-1 


.. 





0\ 




(0 





... 





o\ 





1 


-1 .. 
















... 














1 .. 
















... 


















, Bv — 


















.. 


1 


-1 










... 













.. 









V 


1 


1 ... 


1 


V 



where is an arbitrary real number Ifk^ = the matrices (2.15) define the so-called standard 
boundary conditions. The case j^O corresponds to the "delta potential" of strength (see, 
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e.g., [11]). It is easy to see that the relations 

(kv\ 
kv 



KerAt; = linear span{ 



are valid. Moreover, 



kv 
\kvj 



}, K&vBv = linear span{ 



/ 1 \ 
-1 






1 

-1 





/o\ 




1 

v-v 



} 



RanAj; n RanS^; = RanS^; = linear span{ 



/0\ 





\kv) 



}• 



Thus, independently of the value ofkv, any diagonal unitary Uv leaving both KerA^, and KerBv 
invariant is necessarily of the form Mv = e"^^Ix; with Iv being the deg(i') x deg(o) identity 
matrix. In either case the matrices of this form satisfy (ii) of Theorem I. Therefore, Wt; = Sfor 
all V eV such that \N ^T\'^\. 

Example 2.7. The case A = I and B = provides an example where W(A,B) = U. Indeed, 
since RanA n RanB = {0}, KerA = {0}, and KerS = C"+^'" Theorem I implies that any IX G U 
belongs to \N{A,B). 

We denote the group W from Example 2.6 by Wq. 

Definition 2.8. Wq is the set of all elements U o/ U having the property 



for all V EV. Here the Lv 's are the linear subspaces ofC"'^^"' defined by (1.10). 

We emphasize that the group Wo is independent of the boundary conditions and completely 
determined by the graph S. From Theorem 1 we immediately obtain 

Lemma 2.9. For arbitrary local self-adjoint boundary conditions (A,B) the group Wq is a 
subgroup of\N{A,B). 

The groups Wq and Uo have a common subgroup consisting of elements IX = e"''! with I 
being the unit (n + 2m) x (« + 2m) matrix. Obviously, it is the largest common subgroup, i.e.. 
Wo n Do = S. Thus, the dimension of WqUo is equal to 



(2.16) 



dimWoUo = dimWo + dimUo-l = IVl + pi + 181 - 1. 



3. Magnetic Fluxes 

In this section we will prove Theorem 3 as well as its Corollary 4. The proof uses well-known 
methods of homological algebra. The relevant concepts used below are explained in standard 
textbooks on Algebraic Topology as for example in [33]. 

Let Sint denote the graph obtained from S by removing its external lines. The graph Sim can 
be viewed as a one-dimensional simpUcial complex. 
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We consider the (additive) Abelian groups Ci and Co generated by the 1-simplices a,, / G J 
and the 0-simplices CTp, v eV, respectively, i.e., 



Ci 




For a given internal line labeled by i G J, let t;+ (/) be the terminal vertex which corresponds 
to the endpoint a,- of the interval [0,a/] and V-{i) the initial one which corresponds to the other 
endpoint 0,. 

We define the boundary operator 5i : Ci ^ Co 

Obviously, d\ is a group homomorphism. We, obviously, have 

(3.1) aic = o ^ L L 

(:P+(i)=l' i:V-{i)=V 

for all i; G V(S)- Thus, each element c G Ker is a union of closed loops: c = Y.ieJ^i^i is an 
oriented closed path, not necessarily connected, where the internal Une i is traversed |«/| times 
in the positive direction if n, > and in the negative direction if n, < 0. 

Also we extend the map Ci Co to a chain complex 

O^Z^Co^Ci 

with 

doc = ^ fix) for c = ^ nx)Ox; G Co- 

Obviously we have ^o^i = 0- We call Ker^i the first homology group i/i(S,Z) of the graph S 

and we set //o(S,Z) = Ker do/Ran di to be the zeroth homology group of S. 
With these preparations we give the definition of the flux map. 

Definition 3.1. A map <t> : //i(S,Z) S is called a flux map if it is a group homomorphism, 
i.e., <I>GHom(//i(g,Z),S). 

Here Hom(//i(S,Z),§) denotes the group of all homomorphisms from //i(S,Z) to § with 
unit element 

m:ci-^1gS forall cGi/i(g,Z). 
Observe that //i(S,Z) is a free Abelian group of finite rank. Its rank equals the dimen- 
sion dim//i(S,Z) of //i(S,Z), considered as a Hnear space over Z. Since the flux group 
Hom(//i (S , Z) , §) is die dual of //i (S , Z), it is isomorphic to the torus T'' with d = dimT/i (S , Z) . 
To calculate the dimension d we note that the Euler characteristics ^(Sint) of the simplicial com- 
plex Sint equals \V\ — p|. On the other hand, we have the standard relation £"(9^) = — b\, 
where bo = dim//o(S,Z) and bi = dim//i(S,Z) are the Betti numbers. The complex Sim is 
coimected since S is, so ^jq = 1 • Therefore, 

(3.2) dimHom(//i(g,Z),§) = |a| - \V\ + 1. 

Let IX be an arbitrary diagonal unitary matrix, i.e., U G U. It can be uniquely represented in 
the form 

(3.3) U = diag({ei'P^},e£,{e^'^'"},ea,{e''P" , 

where the sign "— " corresponds to the initial vertex of the internal edge i G J and the sign "+" 
to the terminal vertex. 
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Define the map p : U ^ Hom(7^i (g, Z) , S) via 

p: U^^ui-) 

with 

(3.4) ^u{c) = llexp{ini{(p^ - (pr)}, 
for an arbitrary 

c = ^n,-C7iG//i(g,Z). 

The definition (3.4) is obviously consistent with (1.13) and reflects the additivity of fluxes under 
the union of loops. 

It follows immediately from (3.4) that the map p is a group homomorphism. Recall that U 
and Hom(i?i(g,Z),S) are isomorphic to the tori T"+2'" and with d = dimHi{S,Z), respec- 
tively. Any homomorphism from T"+^"^ to can be represented in the form U i-^ expliMjU 
(component-wise exponentiation) for some {n + 2m) x d matrix M with integer entries (see, e.g., 
[10]). Thus, the kernel of the homomorphism p is a subgroup isomorphic to a torus. 

Lemma 3.2. p is an epimorphism. 

Proof. We have to prove that the map p is surjective. Let Cp,l<p<d with d = AimHi (g , Z) = 
dimHom(i7i(g,Z),S) be an arbitrary basis of generators. It suffices to show that for arbitrarily 
prescribed (e^''' , . . . ,e''^'') G there is a U e U such that p(U) = a>u with given by (3.4) 
which satisfies 

(3.5) ^y^{cp)=e^'', \<p<d. 

All generators Cp are of the form Cp = J^ieJ^^iip)^' with ni{p) G Z. 

Let//i(g,M) be the linear space over R generated by the basis {cpYp^y. We set 

£ = ((pf , . . . , <Pm , • • • , (PmV e M^*". 

Consider the hnear map S : M^*" Hi{S,R) given as 

d 

Obviously, if the equation 

d 

has a solution, then any matrix U defined by <p through equation (3.3) satisfies (3.5). To prove 
that this equation indeed has a solution for arbitrary choice of numbers ^p by the Fredholm 
Alternative it suffices to show that Ker 5^ = {0}. Assume that 

d 

c= Y^apCp^ Ker 5''', G M. 

p=i 

Then 

d 

Y_ OCpfiiip) = for all i G J. 

p=i 

Hence, we obtain 

(d . d 

L <^pni{p))oi = (^pCp = 0, 

iEJ p=l p=l 
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which imphes that c = 0. □ 

Using the First Homomorphism Theorem from Lemma 3.2 it follows that 

(3.6) Hom(//i(g,Z),S)^U/Kerp. 

Obviously, Uq C Ker p. Also one can easily prove that Wq C Ker p. Actually, we have the 
stronger result: 

Lemma 3.3. Ker p = UoWq. 

Proof. First we prove that Wq C Ker p. Let c\,. .. ,Cd be any basis in //i(S,Z). It suffices to 
show that p{\i){cp) = 1 for any 11 G Wq and any basis element Cp. As discussed above, Cp can 
be viewed as a closed loop of length L, i.e. a sequence of triples {v{i),i{i),v'{t)), £= l,...,L 
with t;(1) = v'{L). Let U = diag(e^"i" , . . .,e''P"+^) G Wq with n = |£| and m = p|. Since U G Wq 
all phases <Pj can be labelled by the set of vertices V of the graph S- Obviously, 

P(U)(q) = exp |i Y,{9v'(e) - = exp {ii(pv'(L) - (pv(i})} = 1- 

Thus, we have proved that Wq is a subgroup of Ker p. The inclusion Uq C Ker p is obvious. 

Therefore, UqWo C Ker p. 

Now we calculate the dimension of Ker p. From (3.6) using (3.2) it follows that 

dimKer p = dim U - dimHom(//i (S , Z) , S) 

= |£|+2|j|-|j| + |y|-i 

= |£| + |J| + |V|-1. 

Recall that by (2.16) dimWoUo = |£| + |J| + |V| - 1. Therefore, the groups UqWo and Ker p 
have the same dimensions. Since both groups are isomorphic to tori and UqWo C Ker p, we 
have UqWo = Ker p. □ 

Actually, at this stage the proof of the second claim of Theorem 3 is already completed. 
Indeed, this follows from the relation (3.6) and Lemma 3.3. 

We turn to the case when W = W(A,B) is strictly larger than Wo. Whereas Wo C W(A,B) 
for all local boundary conditions, for non-local boundary conditions it may well happen that 
W(A,B) C Wo (see Example 3.4 below). 

Let d be the canonical quotient map Hom(77i(g;Z),S) Hom(i7i(g;Z),S)/p(W) with 
Ker 6 = p{W). Consider the composition 

p = op. 

Proof of Theorem 3. By construction p is a homomorphism. Recall that Wq C W(A,B) for local 
boundary conditions {A,B). Lemma 3.3 implies that Ker p = UoW(A,B). From Lemma 3.2 it 
follows that Ran p = Hom(//i(g;Z),§)/p(W). Applying the First Homomorphism Theorem 
we obtain the claim. □ 

Proof of Corollary 4. Since Sint is a tree, the group //i(g,Z) is trivial. Thus, the flux group 
Hom(i7i(g,Z),S) is also trivial. Now Theorem 3 implies that U = UoW(A,S). □ 

The requirement of locality of boundary conditions in the statement of Corollary 4 is crucial 
since non-local boundary conditions may act like a loop. This can be illustrated by the following 
example. 
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Fig. 1 . The graph from Example 3.4. It has two vertices, one internal, and two 
external lines. The internal Une has length a, the arrow shows its orientation. 



Example 3.4. Consider the graph depicted in Fig. 1. Consider the magnetic Laplacian with 
the boundary conditions 

AY+B{\j/-\Axir) = 0, 

where 



(3.7) 



A = 



/I 


-1 





"] 

















1 


-1 





B = 




















1 





1 


1 








V 








V 









1 


V 



Herewe use the following ordering: Xj/ = (VAi(0),y/2(0),yA2(a), 1/^3(0))^. The boundary condi- 
tions (3.7) are obviously non-local in the sense of Definition 2. 
It is easy to check that 



KerA = {0}, KerS = hnear span{ 



Ran A = C^, Ran 6 = linear span{ 





1 

/o\ 

1 



















vJ 



/1\ 

-1 





}■ 



Assume there is a gauge transformation G such that the pairs (AU,BU) and {A,B) with 
U = Ug € U define the same Laplacian. By Theorem 1 since U leaves KerA and KerS invariant, 
it has the form 

lX = diag(e'^Sei'^Sei'^,e'^). 

We have 

/O 1 l\ 
11 
110 
\1 1 0/ 

Simple calculations yield [U, A*S]z = 0/or all 



Z&B ^ RanS = hnear span{ 



A*B 









1 












1 
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if and only if = ^ modulo 2k. Therefore the boundary conditions (Alt, SIX) and {A,B) are 
equivalent if and only if 

/ A{x)dx = modulo 2n. 

Jq 

Note that the isotropy group \N{A,B) consists of all 4x4 unitary diagonal matrices which 
are a multiple of the unit matrix. Obviously, W(A,6) is a proper subgroup o/Wq. 
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